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Abstract: We consider curvaton models with potentials that depart slightly from 
the quadratic form. We show that although such a small departure does not modify 
significantly the Gaussian part of the curvature perturbation, it can have a pro- 
nounced effect on the level of non-Gaussianity. We find that unlike in the quadratic 
case, the limit of small non-Gaussianity, I/atl] -C 1 , is quite possible even with small 
curvaton energy density r <^ 1 . Furthermore, non-Gaussianity does not imply any 
strict bounds on r but the bounds depend on the assumptions about the higher order 
terms in the curvaton potential. 
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1. Introduction 

Possible non-Gaussian features of the Cosmic Microwave Background (CMB) tem- 
perature anisotropy can provide important constraints on models of inflation. For 
instance, an observation of a significant non-Gaussianity would effectively rule out 
single-field inflation (for a review of non-Gaussianity, see Usually in the litera- 
ture the non-Gaussianities are characterized by a non-linearity parameter Jnl, which 
is a measure of the non-Gaussian curvature perturbation relative to the Gaussian per- 
turbation. Present WMAP observations yield the limit -58 < /atl < 134 at 95% 
confldence level. With polarization measurements, the Planck Surveyor Mission is 
expected to push the limit down to |/vl| ^ 2.9 0. For the single fleld inflation one 
obtains Jnl which is of the order of the slow-roll pararameters hence a detection 
of non-Gaussianity by Planck would indeed suffice to rule out single-field infiation. 

For the curvaton models p, it has been suggested that a non-observation 
of non-Gaussianity would indicate that the model is ruled out, at least in the case of 
the quadratic potential. In curvaton models the curvature perturbation is generated 
after infiation by the decay of an effectively massless scalar field a different from the 
inflaton. The curvaton energy density remains subdominant until the end of inflation 
so that the density parameter^ r = 4po-/ (4pr + 3po-) ^ 1, where pr is the energy den- 
sity of radiation after inflation. Thereafter the curvaton fleld begins to oscillate and 
behaves effectively like matter. Its relative energy density grows during oscillations 

^Our definition of r is adopted from m 0| and differs from that of e.g pO|. 
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and when it eventually decays, the perturbation it has received during inflation will 
be imprinted on the decay products, the light degrees of freedom. Because the cur- 
vaton is massless, the perturbation is predominantly Gaussian. However, there will 
also be a non-Gaussian contribution that arises because of the curvaton dynamics 
after inflation. By now a well-known result is that in curvaton models with quadratic 



potential the non-linearity parameter can be written as P, ITT 



fNL = l + lr-l. (1.1) 
3 8 3r 

Thus, for low r as required for the subdominance of the curvaton during inflation, 
I/nlI is typically much bigger than 1. However, we wish to point out that this result 
is considerably modified for non-quadratic potentials. Although the curvaton must 
be weakly self-interacting, it is highly likely that there is some departure from the 
quadratic potential. As we will discuss in this paper, while leaving the Gaussian 
perturbation essentially unchanged, a small correction to the quadratic potential can 
have an important effect on the non-Gaussianity parameter. Indeed, we show that if 
one does not insist on a strictly quadratic potential, it is quite possible to have I/nlI 
much less than 1 also in the curvaton scenario. 



2. The curvature perturbation generated in the curvaton model 

We adopt here the non-hnear, the so-called separate universe approach presented 



in |T^, |T^, [T^, [T^, which is valid on large scales. There one considers perturbations 
around the homogeneous and isotropic flat FRW-universe assuming that their spatial 
variation outside horizon is smooth. The evolution of large scales is approximated by 
replacing each quantity by its spatial average inside some smoothing scale k~j^^^^^^ and 
considering these smoothed, locally homogeneous and isotropic regions to evolve like 
separate FRW-universes |l^, The spatial variation of perturbations outside 

the smoothing scale is taken into account by doing flrst order gradient expansion 
leading to different expansion rates in different smoothed regions. Here we briefly 



recapitulate the main features [0, Q of this approach relating to non-Gaussianity 
in curvaton models. 

In the flrst order in gradient expansion, the spatial part of the metric can be 



written as ||T4| 



g,, = a\t)e-'^^''-hij (2.1) 

where is constant, assuming the amplitude of the gravitational waves to be small. 
The curvature perturbation '?/'(t,x) thus deflned can be interpreted as a perturba- 
tion in the scale factor, a(t,x) = a{t)e~'^^^'^\ As it has been shown in [10, 14|, the 



curvature perturbation on uniform energy density hypersurfaces stays constant out- 
side the horizon in the absence of non-adiabatic pressure perturbation; just like its 
counterpart in the usual flrst [O and second order perturbation theories |TB . 



-2- 



The amount of expansion along the worldhne of a comoving observer from a spa- 
tially flat ^jJ = slice at time ti to a generic slice at time t is given by N(t, x) = In^^^^ 
since the expansion in a spatially fiat gauge corresponds to that of the unperturbed 
universe. By choosing the slice at time t to have uniform energy density, the curva- 
ture perturbation on that slice can be written as [|1^, 0] 



at, x) ^ -ij{t, x) = In^^ = N{t, x) - N{t) (2.2) 

where N{t) is the amount of expansion in the background universe. 

Following [|lOl we expand the curvature perturbation ( p. 21) up to second order in 



the Gaussian curvaton perturbations in order to take into account the non-Gaussian 
effects: 

at, x) = d^N{t)6a{t, x) + -dlN{t)6a{t, x)2. (2.3) 

In the limit r -C 1 the curvature perturbation is almost completely generated during 
oscillations of the curvaton field. Assuming sudden decay at t = tdec the amount of 
expansion in the background universe during oscillations is given by flO 



AT/ \ li 1, (Pr)osc /„ -N 

(cTdcc, 0-osc) = -In^^ 2~ = 7 — ^ — (2-4) 

3 2"^^^dec 4 Pdcc - [Pa)dec 

where cTosc is the value of the curvaton at the onset of oscillations. The total curvature 
perturbation is obtained by substituting Eq. (|2.4| ) into Eq. ( ^^ 

where we have denoted 9a-, =' . The initial values for the curvaton field 5(T*(t,x) in 
each smoothed region are set by inflation and the derivatives in Eq. ( |2.3| ) are thus 
taken w.r.t the field value during inflation a^. 

We point out that the nonlinear part in the curvature perturbation, Eq. (p.5| ), 
consists of a square of the linear part (^^5cr*)^ and of an additional dynamical term. 

cr osc 

One could thus expect that the non-Gaussian effects would be more dependent on 
the dynamics than the Gaussian part which is indeed quite reasonable since, loosely 
speaking, the Gaussian part depends on the size of the perturbations while the non- 
Gaussian part depends on the relative size of the perturbations compared to the 
background field. 

From Eq. ( |2.5| ) the curvature perturbation can be written in the form [jl^ 



( = a ~ ^fNiiCg ~ iCg)) where (g is Gaussian and the non-linearity parameter /at^ 
is independent of position. The non-linearity parameter can now be directly read off 
from Eq. (^^) : 



3 



Although this result was obtained in [0, |18|, the dependence of the last term 
'^osc'^osc/o"osc the potential has not been previously examined. In the quadratic 
case cTqsc'^osc/o'osc — but, as we will show shortly, this result may be considerably 
modified even if the deviation from the quadratic potential is small. 

3. Small deviation from the quadratic potential 

The curvaton equation of motion during radiation domination is given by 

a+^^a + V'ia) = (3.1) 

where we have ignored spatial gradients which are small on large scales. We now 
consider the potential of the form 

V{(t) = + Am^-"(T'^. (3.2) 

with A ^ 1. To describe the size of the potential correction at the end of infia- 
tion we introduce a parameter s = 2A((T*/mo-)"~^. The smallness of the correction 
in Eq. ( |3.1| ) requires s <^ 2/n. In the quadratic case the e.o.m ( |3.1| ) is noth- 
ing but a Bessel equation with a general solution ao{t) = AoJ(^i/4^){mo-t) / (mo-tY^^ + 
BQj(^_i/4^){m^t) / {m^ty^^ = AQyi{t) + BQy2{t). To obtain a regular solution at mo-t = 
0, we must set Bq = 0. Furthermore, requiring aolm^jt = 0) = a* we find 

(+\-rr J(l/4)(m^t) 

^o(i) - ^*2V4r(3/4) (m.t)V4 (3.3) 

from which we see that cTosc^osc/cTosc — ^^e quadratic CclSG clS claimed above. 

Let us now make an Ansatz of the form a = ctq + Aai in Eq. ( |3.1| ) . At first order 
in A we obtain the linearized equation of motion 

+ + mlai + m^-^nal-^ = 0. (3.4) 

The solution to the homogeneous equation is already given above; the general solution 
a{t) = A{t)yi{t) + B{t)y2{t) is obtained by the method of variation of parameters. 
The coefficients are solved from the equations 

A'yi + B'y2 = (3.5) 
{A'y[ + B'y'2)ml = -m'-^ nC''-' y^' (3.6) 

where =' and C = o"* 25/4r(3/4) • Since the correction to the quadratic potential 

is small at the end of infiation and even smaller at later times we can assume that the 
beginning of the oscillations takes place at the same time as in the purely quadratic 
case. Thus the era between the end of infiation and the beginning of oscillations 
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corresponds to niat = 0...1. The equations (|3.5|) , ( |3.6|) can now be solved by 
expanding the homogeneous solutions and by a straightforward calculation we find 
that up to order 0{{m^t/2f) 

A{m„t = 1) ^ -n^-^m^-" (^1.0 - O.lOn + 6.9 x 10" V - 

3.6 X 10"V + 1.4 X 10"V) (3.7) 

B{m^t = 1) ^ r2<"^m2-"(^0.82 - 0.095n + 6.8 x 10" V 

-3.6 X 10-V + 1.3 X 10"V) . (3.8) 

Thus the value of the curvaton at the onset of oscillations reads 

cTosc ~ 0.81a* + Anm2-"(T^"-i^(n) (3.9) 

with g{n) = -0.20 + 1.2 x lO^^n - 6.1 x lO^V + 3.1 x lO^V - 2.1 x lO^V. 
With the exponent values < 10 we are considering, g{n) is negative and roughly 
constant, g{n) ~ —0.1. The non-linearity parameter, valid for any potential of the 
type given in Eq. (|3.2|) , is now obtained by substituting Eq. ( ^^ into Eq. ( p.6|) : 



fNL = - + -r--{l + n{n- In - 2) l\ \ J . 3.10 

3 8 3r V (0.81 + 0.50n(n — l)s(7(r2))^/ 

4. The behaviour of I/^vlI in the limit r <C 1 



It is readily seen that the effect of the potential correction in Eq. ( ^.101 ) is most 



significant in the limit of small curvaton energy density r -C 1. In the following we 
are working in this limit if not otherwise stated. Thus we can consider the dominating 
1/r part of /tvl alone and neglect the small contribution from the remaining terms^ 
I + |r. We now examine the effect of the potential correction by keeping s fixed. 
Our solution to the equation of motion (Eq. |3l9|) is constructed in such way that the 
value of the curvaton during infiation, a*, is fixed to be the same as in the quadratic 
case; this is an approximation which however should be justified as we are considering 
only small departures from the quadratic potential. Since we also keep the mass 
unaltered, a constant s means that we are considering the coupling constant as a 
function of the exponent A = \{n). 

As stated above, our perturbative approach to solving the equation of motion 
( |3.1| ) puts limits A -C 1, s -C 2/?7,. Furthermore, the dominance of the Gaussian 
perturbations and the masslessness of the curvaton during infiation also restrict the 



^We will show that the 1/r term may vanish in which case other terms in Eq. (2_^) should 
also be taken into account. The term linear in r is however negligible in the limit r <C 1 and the 
constant part 5/3 does not affect our qualitative conclusions. 
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possible values of a^, and niu. Using Eq. ( p.9|) we find the spectrum related to the 
two-point correlator of the curvature perturbation (i.e. Gaussian spectrum) to be 

n- r:T!^ -f{^ (4.1) 



0.81 + f^(n) 



where {rH^/ (47rcr*)^ is the purely quadratic result; H^, is the Hubble parameter 
during infiation. In the small correction limit the prefactor in Eq. ( |4.1| ) is ~ (9(10^^) 
and hence we conclude that, although the perturbation amplitude is suppressed, the 
Gaussian part is not significantly affected by the correction as the suppression can 
be compensated by a slight increase in the scale of infiation. Thus we may use 
the results for the quadratic parameters which typically imply the restriction 

<ti < coming from the masslessness of the curvaton field and the assumed 
Gaussianity of the curvaton perturbations. This means that the smallness of the 
potential correction, s -C 2/n, requires A <^ 1 when n > 4. Moreover, the infiaton- 
generated curvature perturbation is supposed to be negligible which also requires 
A <^ 1 for non-renormalizable terms, n > 4. 

In Fig. H we show the behaviour of the dominant part of the non-linearity 
parameter as a function of n for two selected values of s. The value of the non- 



Ifnrl 




Figure 1: The qualitative behaviour of the non-linearity parameter with s = 0.2 on the 
left panel and s = 0.05 on the right panel; the values of {/nlI are in units 5/3r = 1. Apart 
from the analytical result of Eq. ( ^.lOD we also show a numerical result obtained directly 
from Eqs. ( |2.6D , ( |3.lD setting t = l/mg- as the beginning of the oscillations. 



linearity parameter depends on the parameters X,a^:,m„, but nevertheless Fig. |T| 
reveals the generic quahtative behaviour of I/tvlI as n is varied. One can clearly see 
that when the exponent of the potential correction is increased, the amount of the 
non- Gaussianity first begins to decrease as compared to the quadratic case. However, 
at large n the value of I/atlI begins to grow rapidly. 

The explanation for the behaviour of {JniI and the physics involved is most 
transparent if we switch to the perturbative point of view. Using perturbation theory 



- 6 - 



one finds [11, IS] that the 1/r part in Eq. ( p.6|) represents the first order contribution 
while the rest is due to second order effects. Thus the first order theory is adequate 
as long as we restrict ourselves in the region where the 1/r term dominates. The 
decrease in the amount of the non-Gaussianity means that the perturbations become 
smaller compared to the background value. In our case this would imply that after 
inflation the perturbations are damped faster than the background value. Indeed, 
we see that this is the case by considering the first order equation of motion for the 
perturbations, 5a + ^5(T + V"{a)5a = 0, with a potential correction of the form 
Am^~"'cr". If the potential is purely quadratic the perturbations and the background 
field apparently obey the same equation. 

The increase of the exponent n diminishes the energy density of the background 
field at the beginning of oscillations since we are keeping s fixed. Furthermore, 
the energy density associated with the perturbations at the end of inflation gets 
bigger. With a large enough n these effects become dominant over the damping of 
the perturbations whence the amount of non-Gaussianity again begins to grow. From 
Fig. |I| we see that the bigger the value of s, the smaller is the value of n at which 
the growth begins; this is of course quite reasonable. 

We should point out here that the increase in I/nlI with a large n typically 
happens when the potential correction becomes significant, s ~ 2/n. At this point 
the perturbative approach breaks down and the result (Eq. (|3.10|) ) can be at most in 
a qualitative agreement with the true behaviour of /atl. The drastic increase in the 
value of I/atlI seen in Fig. |l| is partly due to this effect. In Fig. |I| we also display the 
result of a numerical analysis in which we have, for simplicity, neglected the small 
change in the time tosc corresponding to the beginning of oscillations. The values of 
I/tvlI thus obtained are not significantly different from the perturbative results which 
is to be expected since only the region s < 2/r;, is shown in Fig. The increase in 
I/tvlI levels out, however. 

When the correction s becomes even larger s > 2/n one no longer can ignore the 
effect on the beginning of oscillations. Indeed, the perturbations begin to oscillate 
way before the background field and the oscillation may not initially take place in 
the quadratic part of potential. We do not consider such large corrections in detail 
here but we make some general remarks on the behaviour of f^L justifying the use of 
the perturbative results in the region s ~ 2/n and giving a qualitative understanding 
of the region s > 2/n not covered by our perturbative treatment. Using Eq. ( |3.1j ) 
we find that the condition for a local extremum in fNiin) can approximatively be 
written as 

0"oscO-osc = ^ osc H 2 KoscCnCTosc + (^oscOn(^osc " 2 d^Cr . (4.2) 

To obtain this result we have assumed the potential correction to be of the same 
order of magnitude as the quadratic part. In the limit of small corrections A ^ 1 the 
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expression in parenthesis in Eq. (|4.2| ) vanishes and, since a'^^^ < by Eq. we see 

that in this region /tvl is a monotonously decreasing function of n; this is consistent 
with our perturbative result Eq. ( |3.10| ). However, when the correction becomes larger 
the terms in Eq. (O) involving time derivatives are no longer negligible. Thus, for 
certain values of n there exist solutions to Eq ( |4.2| ) implying that the growth of the 
non-linearity parameter \fNL\ in the region s > 2/n eventually ceases whereafter 
I/atlI begins to oscillate. We do not examine the non-Gaussianity in this region 
more closely but we point out that I/atl] might become large enough to exclude 
some classes of potentials already with the present WMAP limits 0. 



5. Restrictions on the potential correction 



So far we have considered the potential corrections only from a technical point of 
view. There are, however, some physical motivations for choosing non-quadratic po- 
tential. Small corrections of the type n = 2 + t typically represent the effects coming 
from one-loop corrections due to light degrees of freedom the curvaton couples to. 
As we have seen above, these tend to decrease the amount of non-Gaussianity. Also, 
it is conceivable that the curvaton is self-interacting. The n = 4 term, in particular, 
would be interesting since it implies less fine-tuning to satisfy the smallness condition 
of the correction s ^ 2/n; for n = 4 one would not have to require A <^ 1 as for 
the higher order cases such as might arise, for example, in models [^,21] where the 
curvaton field is considered to be one of the flat directions of the Minimally Super- 
symmetric Standard Model. This is also seen in Fig. ^ where we show /tvl as a 
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Figure 2: The non-linearity parameter \fNL\ in units 5/3r = 1 as a contourplot with 
n = 4 on the left panel and n = 6 on the right panel. The scale on y-axis is logarithmic. 
The values of I/atlI are evaluated only in the perturbative region s < 2/n and the non- 
perturbative region is printed in white. 
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contourplot in [a jm^, A) space for n = 4 and = 6. We note that especially in the 
n = 4 case there is a significant region in the parameter space in which the poten- 
tial correction is small s < 2/n but the value of /tvl is highly suppressed from the 
quadratic case. For non-renormalizable terms (n > 4) the requirement of negligible 
inflaton-generated curvature perturbations sets an upper limit on A (e.g. for n=6 
A < 10^^°), but the region in the parameter space with small values of I/atlI is still 
considerable. In other words, it is quite possible to obtain I/atlI -C 1 even in the 
curvaton models by adding a small self-interaction term to the quadratic potential. 

WMAP yields an upper limit [|| \ jNL\ ^ 100, which in the case of a quadratic 
curvaton potential implies that r > 0.02 as it can be seen from Eq. (|1 . 1| ) . For the 
non-quadratic case, the limit is greatly modified due to the decrease in f^L. From Eq. 
( P.10|) the part of the parameter space (s, r) compatible with present observations is 
given by 

n{n - l){n - 2)g{n){0Als + 0.25s'^ng{n)) 



1 

r > — 
60 



(5 I) 

{0.81 + 0.50n{n-l)sg{n)y ' ^'^ 

The allowed region is represented in Fig. 13 for a choice of parameter values. It is 




noteworthy that the limits on r are strongly dependent on the size s and form n 
of the potential correction. In the non-perturbative region not shown in Fig. |^ we 
expect the increase in the lower limit on r to level out as a result of the behaviour 
of /nl described above. 

6. Conclusions 

In this paper we have examined how a small departure from the quadratic curvaton 
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potential will affect the produced level of non-Gaussianity of the primordial curvature 
perturbation. For non-quadratic potentials there are two competing, opposite effects 
that contribute to the net non-Gaussianity. First, unlike in the quadratic case, the 
perturbations will be damped faster than the background value. This tends to reduce 
the value of the non-linearity parameter |/iVL|. Second, the increase of the exponent 
n diminishes the energy density associated with the background field at the beginning 
of oscillations. With for a steep enough potential this effect becomes dominant and 
compensates for the damping of the perturbations. As a consequence, the value of 
I/tvlI will increase as compared with the purely quadratic case. 

The net outcome is that although a small departure from a quadratic curvaton 
potential does not modify significantly the Gaussian part of the perturbation, it 
can have a pronounced effect on the level of non-Gaussianity. In particular, the 
limit I/atlI -C 1 is allowed even with small curvaton energy densities r <^ 1. This 



is in sharp contrast with the quadratic result p, |Tl| and shows that the curvaton 
models are not ruled out by a possible non-detection of non-Gaussianity as it has 
been suggested in e.g. [§. Furthermore, unlike in the quadratic case, there are no 
strict limits that could be placed on the curvaton energy density parameter r. By 
adding a small correction to the potential, in practice one can enable arbitrarily small 
r with suitably chosen parameter values. This is an interesting result in the sense 
that it implies that one can not use present observations to fix the lower limit for 
the energy scale of an approximately quadratic curvaton potential without making 
further assumptions about the higher order terms. 
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